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ABSTRACT 



The 3D magnetic topology of a solar active region (NOAA 10956) was recon- 
structed using a linear force-free field extrapolation constrained using the twin 
perspectives of STEREO. A set of coronal field configurations was initially gen- 
erated from extrapolations of the photospheric magnetic field observed by the 
Michelson Doppler Imager (MDI) on SOHO. Using an EUV intensity-based cost 
function, the extrapolated field lines that were most consistent with 171 A pass- 
band images from the Extreme UltraViolet Imager (EUVI) on STEREO were 
identified. This facilitated quantitative constraints to be placed on the twist (a) 
of the extrapolated field lines, where V x B = aB. Using the constrained values 
of a, the evolution in time of twist, connectivity, and magnetic energy were then 
studied. A flux emergence event was found to result in significant changes in the 
magnetic topology and total magnetic energy of the region. 

Subject headings: Sun: Magnetic Fields - Sun: Flares - Sun: Corona 
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Introduction 



Extreme solar events such as solar flares and coronal mass ejections (CMEs) originate 
in the corona of active regions. As the corona is a low-/? plas ma, the mag netic field 



rather than the gas pressure governs its structure and evolution ([Gary. 



200 lj ). This enables 



> of the photosp 


reric 


(Longcope et al. 


2009) 



reric field to be used to study the detailed 



20091 ) . Of particular interest to space weather 



applications is how the three-dimensional (3D) field topology of active regions governs the 
production of extreme solar events. Unfortunately, previous attempts to reconstruct the 3D 
topology of active region fields have been limited by observational issues, such as incomplete 



measurements of the photosph eric vector fie 
extrapolation algorithms (e.g.. iDeRosa et al. 



d, or theoretical issues such as ill-constrained 



2009 



and references within). With the launch 



in December 2006 of the Solar TErrestrial RElations Observatory (STEREO) satellites, the 
observational limitations on reconstructing the 3D topology of active regions has largely 
been overcome, thus opening a new era in the study of the ever-changing solar corona. 

Observations of the coronal magnetic field are difficult to obtain due to the tenuous 
nature of the corona and the absence of strong magnetically sensitive emission lines. 
Estimates of the strength of the coronal magnetic field are possible using radio observations, 
however, this remains a difficult problem requiring simultaneous knowled ge of both the 



temperature and density of the solar atmosphere (jBrosius fc White 



20061). In order to 



overcome these issues, modellers have developed a yariety of theoretical techniqu es to 



examine the topology of the coronal magnetic field (IGary 



1989 



Mikicetal 



19901 ) . These 



methods require observations of the photospheric or chromospheric magnetic field as 
boundary conditions for the extrapolations. Well constrained extrapolations are norma 



produced using photospheric vector magnetograms, although 



Metcalf et al. 



(1995 



20021 ) 



ly 



have developed methods that make use of chromospheric vector magnetograms. As the 
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chromosphere is in a force-free state, chromospheric vector magnetograms are ideally suited 
as boundary conditions in force-free extrapolation methods. However, these observations 
have a limited field-of-view and spatial and temporal resolution. As such, photospheric 
magnetograms are commonly used as boundary conditions for magnetic field extrapolations. 

The extrapolation of the photospheric field into the corona requires that a number of 
assumptions be made. Firstly, it is assumed that the coronal magnetic field remains stable 
for the duration of the observation. Therefore, a static magnetic field model is normally 
used. In addition, it is reasonable to assume that the field is in a force-free configuration. 
This results in a vanishing Lorentz force, j x B = 0, where j is the current density, B is the 
magnetic field. The force-free condition requires that the current density remains parallel 
to the magnetic field, /^oj = aB, and therefore V x B = «B, where a is the amount of twist 
in the field and is some times referred to as the torsion function. In the case of a = a(r), 
the governing equations becomes V x B = a(r)B and B • V«(r) = 0. This non-linear 
force free (NLFF) field is the generic solution to the force-free field condition. Complete 
knowledge of the photospheric magnetic field allows for the inversion of the preceding 
equations and the calculation of a at each point on the photospher e. As such, vect or 



2006|). 



magnetograms are required to calculate the NLFF field in the corona (jSchrijver et al. 
While NLFF methods offer the possibility of increased accuracy, they are computationally 
intensive. Additionally, NLFF extrapolation suffer from several other sources of error, such 
as the limited field-of-view, non-force-free nature of the photospheric field, the intrinsic 
180-degree amb iguity, and low sign al to noise of the transverse component of vector 



magnetograms (IDeRosa et al. 



20091 ). Forcing a = results in the potential or current-free 
solution, which is a first order approximation to the coronal magnetic topology. The 
potential field solution is the first of two subclasses of solutions to the force-free field. Given 
the reduced complexity of the potential solution, the only boundary condition necessary 
to form a solvable set of equation are measurements of the line-of-sight (LOS) magnetic 
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field (lGaryl I1989T ) . The second and final subclass of solution to the force- free field is the 
linear force-free (LFF) field, where a = constant and V x B = qB. Assuming that a is 
constant for all field lines allows for the formation of a closed set of equations wit h only 



knowledge of the LOS magnetic field, provided a is predefined ( lAlissandrakis 



19811 ). 



A major difficulty in studying the extrapolate d field line distribut i on is deciding on 



what is the most appropriate 



recently 



Burnette et al 



way to constrain a. 



Leka & Skumanichl (119991 ) and more 



(j2004f ) compared three methods and applied them to a single active 
region. The first method considered was based on the best fit value of a from a LFF 
extrapolation of a longitudinal magn etogram that best m atched the horizontal field of the 



corresponding vector magnetogram (jPevtsov et al. 



19951 ). Sec ondly, the value of a was 



chosen from the mean of the distribution of local a z = a(x,y) (jPevtsov et al.lll994l ). In the 



final method, the value of a = J z /B z was inferred fro m the slope of a 



Burnette et al. 



east-squares linear 



(12004 conducted a 



fit to the distribution of local J g (x,y) versus B z (x,y). 
detailed study of these three methods, concluding that they give statistically consistent 
values of a. As the methods described above require the availability of high quality vector 
magnetograms, they are not always an available option for constraining the values of a. 

EUV and X-ray images of emitting coronal structures, mainly loops, can also be 
used to constrain the choice of a in magnetic field extrapolations. Initially, loop tracing 
methods were used to provide a first guess at the structure of coronal loops and throug h 



the application of numerical methods the 3D structure found (lAschwanden et al. 



19991 ) 



Other authors have developed methods that minimize the separation between the field lines 



calcu la ted from observations and e xtrapolations a s a con straint on the choice of a ( 



2003). 



Wiegelmann et al. 



(J2005J) 



Carcedo et al. 



(12003 ). and 



Aschwanden et al. 



im et al. 



(12008a 



developed related techniques that examine the EUV emission of coronal 
retrieve their two dimensional structure using a cost function method. 



loop structures an d 



Wiegelmann et al 
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( 120051 ) developed a cost function grading method to evaluate the agreement of calculated 
field lines and the observed emission of coronal loops. This enabled them to gauge 
the optimum value of a for a small number of LFF extrapolations. The cost function 
method was found to allow for the selection of the optimum value of a used in the LFF 
extrapolations, however, the number of possible values of a e xamined was insuffi c ient to 



allow for a detailed study of the active region. Most recently, lAschwanden et al. 



(2008b 



20091 ) expanded on the methods discussed above with the use of the STEREO spacecraft 



to study the structure and oscillations of coronal loops. 

Once the value of the magnetic twist (a) in a region has been well constrained, an 
accurate estimate of the magnetic energy or more importantly, the free magnetic energy, can 
be estimated. The magnetic energy budget of active regions can be calculated using volume - 



integ ral and surface-integral methods such as the magnetic Virial Theorem (jEmslie et al. 
2004 ). The free-energy is calculated as the difference between the magnetic energy of the 
NLFF field and that of a lower bound such as the potential (current-free) field. These 
methods have shown the free energy i n active region's to be in the range of 10 32 — 10 33 ergs 



depending on the size of the region (IGary 



1989 



1995 



Regnier fc Priest 



2007 



Metcalf et al. 



20051 ). The magnetic Virial Theorem assumes the measured field is a force-free field 
and is only possible if the measured data is force-free consistent. Additionally, knowing 
the topology of the active region at instances during its evolution, changes in the regions 
connectivity can be examined. Studies of the changing connectivity of active regions using 
a combination of act ive region segm e ntation techniques and magnetic field extrapolatio ns 



Longcope et al. 



have been achieved (|Longcope 



2001 



Longcope &: Beveridge 



2007 



Longcope et al.l 120091 ) 



(120091 ) concluded that comparison of the changing connectivity in active 



region over time could provide insight into energetics and reconnection in the coronal fields. 



The methods described in this paper allow for the 3D structure of the magnetic field 



-7- 



of an active region to be recovered and a detail ed analysis of its evo lution to be performed 
using SOHO Michelson Doppler Imager (MDI, 



Scherrer et al 



Howard et al 



1995) LO S magnetograms 



200Q ) images. The 



and STEREO Extreme Ultra Violet Imager (EUVI, 
observation and the transformation of co-ordinates between STEREO and SOHO reference 
frames are outlined in Section while details of the extrapolation method and the cost 
function are then given in Section [3J Our results, relating to twist, connectivity and 
magnetic energy are presented in Section HI and our conclusions are given in Section [5j 



Observations 



NOAA 10956 produced a GOES class B9 flare on 2007 May 19. Associate d with this 
flare was a CME and a globally propagating disturbance within the solar corona (ILong et al. 
20081 ). Extrapolations of the region's coronal magnetic field were constrained at the lower 
boundary with LOS magnetic field measurements obtained by MDI on SOHO. MDI images 
the Su n on a 1024 x 1024 pixel 2 CCD camera through a series of increasingly narrow 



filters ( Scherrer et al. 



19951 ). A pair of tunable Michelson interferometers enable MDI to 
record filtergrams with a FWHM bandwidth of 94 mA. 96-minute cadence magnetograms 
of the solar disk with a pixel size of ~ 2" are used in this work. 



Images from EUVI 



Investigation (SECCHI, 



part of the Solar 



Howard et al. 



Carth Connection Coronal and Heliospheric 
20001 ) suite of instruments on board both STEREO 
(A) head and (B)ehind, were used in order to constrain the choice of a within the LFF 
extrapolation of the corona. In particular, 171 A images from EUVI were used as this 
passband has few er contributions from diffuse higher temperature emission lines compared 



to 195 A images (jPhillips et al.l 



20051 ) . This results in more defined loops which are easier 
to compare by the methods described in Section [3j Each STEREO image was co-aligned 
with a common SOHO MDI image using a transformation method based on the work 
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of 



Aschwanden et al 



( I2008bl ). so as to allow the transformation of points between the 



STEREO A B and SOHO reference frames. 



3. Methods 

Given the limited reliability of vector magnetograms and the errors associated with 
their use in NLFF methods, a LFF methods is used in this study to recover the coronal field 
as accurately as possible. Assuming a force-free field, j x B = 0, a LFF field is governed by: 

VxB = «B, (1) 

where a = constant. For LFF field extrapolations, a is a predefined free parameter. Small 
changes in a result in large changes to the extrapolated field. Comparing extrapolated 
fieldlines with EUV observations provides a method of restricting a and improving the 
accuracy of our results. 

It is possible to systematically determine the value of a that provides the best possible 
comparison with the EUV observations using the following steps: 

1. Compute the linear force- free field. 

2. Calculate a large number of possible magnetic field lines. 

3. Project field lines onto EUV observation. 

4. Calculate the cost function along each field line. 

5. Repeat for each additional value of a. 

6. Compare the cost of the best matched field lines for each a. 




Fig. 1. — Field lines overlaid on emission with a Gaussian profile. Top: Bright Gaussian 
ridge, with three overlaid field lines 'A', 'B', and 'C. Field line 'A' crosses the ridge, 'B' rest 
along the ridge, and 'C oscillates along the ridge. Bottom: Intensity as measured along each 
field line. 
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In order to do this accurately it was first important to define what a good fieldline is. 
Figure [T] shows the path and resulting intensity for three different ideal field lines. Field 
line 'A' crosses the bright emission only once and has a small spike in emission along its 
length. Field line 'B' rests along the field line and has a smooth bright intensity profile. 
Field line £ C oscillates on and off the bright emission and the resulting profile has a series 



of spikes in the intensity. With this in mind a number of methods have b een developed 



to 



systematically extract loop pr ofiles from EUV image of the solar corona (jCarcedo et al 



2003 



Wiegelmann et al 



20051 ) . Including this information in the definition of a cost function 



used to constrain the choice of a removes the intermediate step of extracting coronal loo p 



Wiegelmann et al. 



fl2005h 



profiles. For the purposes of this research, modifications to the 
method were studied to find the most computationally efficient and best performing cost 
function. A small sample of the various cost functions studied are: 



C w = jVI(l)dl/ (l(j I(l)dlf 
C EW = jVI(l)dl/ (l j I(l)dl 
C B = l/fl(l)dl, 

where 1(1) is the intensity of emission along a loop and I is the 
are the standard and equal weight cost functions as defined by 



loop length, Cw a n d C_r w 



Wiegelmann et al. 



(120051 ). 



and Cb is a modified cost function developed to select the brightest field lines. 



Using the LFF method from 



Alissandrakisl ( Il98ll ). field lines were extrapolated starting 



from a user-defined foot-point. Foot-points were selected from EUV observations of the 
coronal loops in question. Field lines were then transformed to the desired perspective 
of each STEREO spacecraft, and the emission along each field line and associated cost 
calculated. Figure [2] ( Top) show the ability of the three cost functions to recover the 
geometry of the active region loop. As can be seen, the Cew an d Cw cost functions perform 
best, both however include unrealistic field lines. These are due to the infinite number 
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Fig. 2. — Field lines that best represent the observed emission as calculated by the selection 
algorithm for a = 6.25 x 1CT 3 Mm^ 1 . (a) Cw, (b) Cew an d (c) Cb cost functions without 
field line path restrictions, {d, e, f) same with field line path restriction as outlined with 
overlaid circles. 

of possible field line paths leaving any sub-region of the active region. Given the tenuous 
nature of the coronal loops we wish to study, field lines that cross the core of the region and 
other high emission locations will have unacceptable paths and relatively low costs. 

In order to overcome this, fieldlines were forced to emerge from user defined foot-points 
and pass through a number of guide circles. These user defined guide-circle domains are 
highlighted in Figure |2] {Bottom). Figure [2] [Bottom) shows the ability of the three cost 
functions to recover the geometry of the active region loop with the inclusion of these 
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Fig. 3. — (Top) 15 field lines with the smallest cost as selected using the Cew cost function 
with STEREO A EUVI images of NOAA 10956. (Bottom) Same but for STEREO B images, 
for a = 8.7 x 1(T 3 Mm" 1 . 



guiding circles. This can be seen to greatly increase the ability of each cost functions to 
recover the coronal geometry. The Cew cost function was used in this work to select the 



Fig. 4. — 3D view of field lines shown in Figure [3j 



value of a for the LFF extrapolations due to its ability to select more valid field lines under 
both conditions. The user-defined foot-points and guide circles stop the algorithm from 
selecting loops and other bright features outside the region of interest. Large guide circles, 
with radii less that a typical loop half-length, were chosen to ensure that difference user 
selections would result in similar features being used to constrain the extrapolations. It 
may, in the future, be possible to automate the selection of coronal loops based on EUV 



and/or X-ray emission, but this will 
coronal loop detection algorithm (e.g 



require the application o 



Aschwanden et al. 



2008ah . 



an accurate and robust 
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Using the Cew cost function it is possible to compare the cost function for various 
values of a and determine the value that most closely resembles the observed field. More 
than 5,000 field lines are plotted for each value of a and the cost along each is calculated. 
In order to enhance the sensitivity of the method to changes in the values of a only the 
total cost of the 15 best field lines are compared. As shown in Figure EJ this is a sufficient 
number of field lines to accurately recover the topology of the region. The contrast of the 
cost for each value of a was increased further by normalizing the returned values and setting 
the cost to 1 for values of a where no field lines pass through the guiding path and/or end 
at the user-defined foot-point. 

Figure [3] shows the best field lines as selected by the method for STEREO EUVI 
images from both the A and B spacecrafts. The field lines returned by the method are not 
the same for each spacecraft. However, the general geometry defined by the field lines is 
very similar. Figure H] shows the field lines from both perspectives plotted together in 3D. 

An estimate of the energy of the LFF field in excess of the potential field can be 
obtained from: 

Emjf = E LFF — Ep ot , (2) 

where Elff is the magnetic energy of the LFF field extrapolation, 

Elff = J (B 2 LFF /8n)dV, (3) 

and E po t is the energy of the potential field extrapolation, 

E Pot = J (B Pot /8n)dV, (4) 

where dV is a volume within the computational domain, and Ep ot and Elff are the 
magnetic energy of the potential field and LFF field respectively. 

As a precautionary warning, it should be noted that LFF extrapolations have a 
maximum value of a that can be studied before oscillatory signals from the complex 
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solution of the governing set of equation begin to dominate the systems (iGaryl 119891 ). We 
have investigated this effect for the data set studied here and found that for a values up to 
0.03 Mm^ 1 the oscillatory signal from the complex solutions is negligible up to a height of 
100 Mm, the maximum height of loops studied here. 



4. Results 

Using the method described in Section [31 an initial investigation of the cost function 
behavior for a range of a values from —0.3 Mm -1 to 0.3 Mm -1 had established that a 
minimum in the cost function existed and had a positive value of a. To verify this a more 
detailed analysis of the cost function was conducted for twenty three values of a within the 
range of —0.11 Mm -1 to 0.28 Mm -1 . Figure [5] shows the evolution of the cost function for 
these values of a at the times studied. At all times there is a strong correlation between 
both the STEREO A and B cost functions. At 09:35 UT there is a distinct minimum in the 
cost at 5 x 10 -3 Mm -1 . The cost for the STEREO B observations was seen to decrease for 
larger values of a, however as the STEREO A cost does not follow this trend it is assumed 
false. At 11:11 UT the minimum has increase to around 1 x 10 -2 Mm -1 , the relatively high 
cost returned for a = 5 x 10 -3 Mm -1 reinforces this increase. From 12:47 UT onwards the 
minimum cost in a or twist of the coronal loop analyzed is seen to decrease to levels seen 
at 09:35 UT. 

The evolution of magnetic flux within each sub-region was examined in order to more 
fully understand these changes in a. Each magnetogram image of NOAA 10956 was 
thresholded at the ±100G level to identify the main regions of magnetic flux in the global 
structure. Figure E] shows the main positive and negative sub-regions within the active 
region, ranked by area from largest to smallest. From the corresponding STEREO /EUVI 
image it is clear that the large loop structure studied here is associated with field lines 
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connecting regions 1~ and 3 + . Table [T] summaries the changing magnetic flux within the 
region (some negative regions are excluded due to changes in ranking during the time 
studied). 

As can be seen there was a significant changes in the magnetic flux in the region during 
the time analyzed. From 09:35 UT to 11:11 UT, the magnetic flux in region 1 + increase by 
0.6 x 10 20 Mx, while in regions 1 and 2 + the magnetic flux decreased by 1.3 x 10 20 Mx 
and 1.1 x 10 20 Mx respectively. During this time the other sub-regions saw little or no 
changes in their magnetic flux. From 11:11 UT to 12:47 UT, the magnetic flux in regions 
1 + and 2 + increased by 0.6 x 10 20 Mx and 0.3 x 10 20 Mx respectively. While the magnetic 
flux in region 3 + decreased by 0.3 x 10 20 Mx. From 12:47 UT to 14:27 UT, the magnetic 
flux in regions 1+, 1", and 3+ increased by 0.8 x 10 20 Mx, 0.4 x 10 20 Mx, and 1.8 x 10 20 Mx 
respectively. 

Knowing the values of a, it is then possible to study changes in the connectivity within 
the region. As region 1~ is the common foot-point for most of the loops in the active region 
it was selected as the source for all the extrapolated field lines used to examine the region's 
changing connectivity. Table [2] outlines the changing connectivity for the time studied. 
From 09:35 UT to 11:11 UT the major changes in connectivity are associated with the 

Table 1. Evolution of the magnetic flux in the major sub- regions contained within NOAA 

10956. Results are reported in 10 20 Mx. 



Region 09:35 UT 11:11 UT 12:47 UT 14:27 UT 



1- 


77.8 


76.5 


76.5 


76.9 


1+ 


32.5 


33.1 


33.7 


34.5 


2+ 


12.5 


11.4 


11.7 


11.9 


3+ 


7.2 


7.1 


6.8 


8.6 


4+ 


3.3 


3.4 


3.7 


2.1 
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transfer of flux from a closed to open configuration and a 1.9% increase in flux connecting 
region 1~ to region 4 + . From 11:11 UT to 12:47 UT, there is a decrease in the amount 
of open flux and an increase in the connectivity from region 1~ to regions 1 + , 2 + , and 
3 + by 4.8%, 0.8%, and 2.3% respectively. Additionally, the amount of flux joining region 
1~ to region 4 + decreases by 3%. From 12:47 UT to 14:27 UT the amount of open flux 
continues to decrease and there is no flux left connecting regions 1~ to 4 + . The amount of 
flux connecting region 1~ to region 1 + increases by an additional 9%. There are marginal 
changes in the flux connecting 1~ to 2 + and 3 + . 

Following from the investigation of the changing magnetic flux within the region and 
the significant changes in connectivity, an investigation was undertaken into the energy 
stored in the magnetic field. Table [3] shows the calculated magnetic energy for two volumes: 
one centered on the loop structure connecting regions 1~ and 3 + (see boxed region outlined 
in Figure [6]) and the other being the full computational volume. An estimate of the error in 
the energy calculations is provided by calculating the energy of the LFF field for a m i n — A a 
where A a = O.lTu and a m i n is the calculated a value of the loop observed. From this it is 
clear that there is a significant buildup in energy within both volumes from 09:35 UT to 
11:11 UT. At 12:47 UT, the difference in magnetic energy between the potential and LFF 
field has significantly decreased and is only marginally higher than that at 09:35 UT, while 
at 14:27 UT this difference in magnetic energy is lower than at 09:35 UT. 

5. Conclusions 

NOAA 10956 rotated onto the solar disk on 2007 May 15 and on 2007 May 19 the 
region produced a B9 class flare, starting at 12:34 UT, peaking at 13:02 UT, and ending at 



x The step size of a during the cost function analysis is A a = 0.17 
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Table 2. Connectivity matrix for field lines starting at region 1 , with threshold of 
±50 G. All results are a percentage of total flux leaving region 1~. 



1" 


09:35 UT 


11:11 UT 


12:47 UT 


14:27 UT 


Opcn 


13.9 


20.6 


15.9 


12.9 


Closed 


31.9 


24.3 


23.8 


23.6 


1+ 


24.5 


25.1 


29.8 


38.9 


2+ 


11.1 


10.8 


11.6 


11.2 


3+ 


9.8 


9.7 


12.0 


12.3 


4+ 


7.4 


9.3 


6.2 


0.00 



Table 3. Changes in the magnetic energy contained within different volumes in NOAA 
10956. The small region is a box centered on the the loop structure connecting regions 1~ 

and 3 + . The full volume in the whole computation domain. All values for energy are 
reports in units of 10 32 ergs.. Associated errors in the calculation of the LFF energy and 
the energy difference are on average ±0.01 (±0.03), and a maximum of ±0.03 (±0.05) at 

11:11 UT for the small region (full volume). 



Size 


Parameter 


09:35 UT 


11:11 UT 


12:47 UT 


14:27 UT 


Small Region 


Pot 


4.46 


4.44 


4.44 


4.48 




LFF 


4.48 


4.50 


4.46 


4.49 




Difference 


0.013 


0.061 


0.021 


0.01 


Full Volume 


Pot 


6.51 


6.38 


6.33 


6.23 




LFF 


6.53 


6.49 


6.37 


6.24 




Difference 


0.021 


0.109 


0.037 


0.018 
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13:19 UT. The flare location suggests that the source of the flare was the loop structure 
connecting regions 1~ to 2 + in our analysis. According to the SOHO /LASCO CME 
Catalogue^ , the region had an associated CME which first appeared at 13:24 UT in the C2 
field of view. In addition, a disappearing filament was observed eru pting from the lower 



right section of the regions beginning at 12:31 UT ( lLong et al. 



20081). 



Our analysis has highlighted a number of physical changes within the magnetic 
structure of the active region around this time. These are: 

• Magnetic flux cancellation and emergence. 

• Increase in the amount of a or twist along the observed loop. 

• Changes in the connectivity between sub-regions. 

• Significant changes in the magnetic energy stored in the region. 

The increase in the amount of twist or a along the observer loop from 09:35 UT to 
11:11 UT, was accompanied by an increase in the magnetic flux within sub-region 1 + , 
and decreases in the magnetic flux within sub-regions 1~ and 2 + and an increase in the 
amount of magnetic energy stored within the volume. The increase in magnetic energy can 
be explained by the increase in a used to calculate it. However, as a is obtained from a 
comparison to EUV observations this increase is taken as real. Following the increase in 
a at 11:11 UT, the amount of twist along the observed loop decreases to 09:35 UT levels 
by 14:27 UT. During this time period the magnetic flux in sub-regions 1 + , 1~, 2 + , and 3 + 
increases by 0.4 — 1.5 x 10 20 Mx and there is a 1.1 x 10 20 Mx decrease in the magnetic 
flux in region 4 + . Additionally, the magnetic energy decreases to 09:35 UT levels over the 
same time period. Changes in the connectivity from region 1~ can be associated with the 



http : // cdaw . gsf c . nasa . gov/CME_list 
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changes in a prior to 11:11 UT. However, by 12:47 UT the amount of flux connecting to 
region 1~ to regions 1 + and 3 + rose significantly. 

From 09:35 UT to 11:11 UT the difference in magnetic energies increased from 
1.3 x 10 30 to 6.1 x 10 30 ergs in the sub-volume centered on the observed loop and from 
2.1 x 10 30 to 10.9 x 10 30 ergs in the full computational domain. This result is particularly 
important as the energy calculation at 11:11 UT is significantly different from the potential 
case, as indicated by the error estimate. It should be noted that smaller error bounds may 
be possible provided that a more detailed analysis of the cost function (smaller A Q ) is 
performed and that a clear minimum is the cost function exists. By 12:47 UT this difference 
in magnetic energy in the system had dropped by a factor of three and by 14:27 UT had 
decreased below the difference at 09:35 UT. The build-up and subsequent decrease in 
magnetic energy in the system is correlated with the mechanisms thought to be responsible 



for the release of the B9 flare at 12:30 UT. The method is unable to calcu 



ated the free 



19781 ). However, the 



energy of the system due to the limitations of LFF methods (ISeehaferl 
reported changes in the difference in magnetic energy between the potential and LFF fields 
provide a estimate for changes in the free energy contained in the region during this time. 

The above investigation of magnetic properties of NOAA 10956 has shown that from 
09:35 UT there was a significant amount of flux cancellation within the region, accompanied 
by an increase in the amount of twist and free energy along the observed loop. Following 
the onset of the flare at 12:34 UT, the amount of twist and magnetic energy in the loops 
decreased. This was accompanied by the decrease in open flux leaving region 1~ and 
increase in the magnetic flux in the region. It is suggested that the sudden decrease 
and subsequent increase in magnetic flux within the active region and the changes in 
twist and magnetic energy caused small changes in the observed connectivity between 
sub-regions, Table [2j The emergence of magnetic flux within the region and resulting 
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increase in magnetic energy could have resulted in the region passing a critical threshold, 
passed which it was in an unstable configuration. Under the theory of self-organized 
criticality this would have resulted in a s eries of events returning the region to a stable 



configuration (IVlahos fc Georgoulis 



20041 ) . Whether the mechanism responsible for passing 
this critical threshold is the unbalanced emergence of flux within the region or the resulting 
non-potentiality of the flux ropes remains unknown. The decrease in both a and magnetic 
energy and increase in connectivity of loop connecting region V~ to region 3 + suggest that 
the loop had room to expand and relax into the neighboring sub-region following the flare 
and eruption in the neighboring sub-region. 

The method described in this paper allows for the detailed analysis of the amount of 
current or a in coronal loop structures and with the added use of the t win perspectives of 



STER EO EUVI images, are an 



( 120051 ) and 



Carcedo et al 



expansion of the methods developed by 



Wiegelmann et al. 



( 120031 ). These previous studies involved the use of cost functions 
as a mechanism for constraining the free parameters in extrapolations of the solar corona 
and were restricted to the investigation of one or two properties and singular instances in 
time. More advanced NLFF method have been used to investigate the free magnetic energy 



of active r egions (jRegnier fc Priest 



2007; 



evolution (IThalmann fc Wiegelmann 



Schriiver et al. 



2008 



Thalmann et al. 



2008) and in some instance its 



20081 ) . However these methods 



have been restricted by the limited time cadence of vector magnetogram data and provide 
little detail into the physical changes in the structure of active regions surrounding the 
buildup and release of magnetic energy. While the method developed here are unable to 
accurately calculate the free energy of the active region, changes in the energy difference 
can be used as a proxy to changes in the free energy. 

The method has been shown to detect changes in the amount of twist (and hence 
current) within coronal loop structures using the twin perspective of STEREO. Coupling 
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the temporal analysis of the structure of NOAA 10956 with an investigation of the magnetic 
flux within each sub-region has provided great insight into the evolution and driving forces 
within the active region. Future work is needed improve the method and remove the need 
for user input. Expanding the method for the analysis of multiple loop structures in active 
regions would greatly increase the diagnostic power of the algorithm. 

We are grateful to the referee for their helpful comments. P. A. Conlon is an IRCSET 
Government of Ireland Scholar funded under the Irish National Development Plan. The 
authors would like to thank Dr. R.T.J. McAteer and Dr. D.S. Bloomfield for their helpful 
discussion and thoughts. The LFF method used in this study is an adaptation of code 
originally developed by Dr. Thomas R. Metcalf. 
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Fig. 5. — Top to bottom: Cost for a range of a values for MDI magnetograms obtained at 
09:35 UT, 11:11 UT, 12:47 UT, and 14:27 UT. Solid and dashed lines represent comparisons 
with STEREO A and B EUVI observations respectively. Error bars in top panel represent the 
associate error in all the values of a due to the step-size used in the calculations. Minimum 
in the cost function and the associated error are highlighted by the vertical dotted lines 
and dashed-dotted lines respectively. The active region produced a B9 flare that began at 
~12:34 UT, peaked at ~13:02 UT, and ended at -13:19 UT. 




Fig. 6. — Sub-regions of NOA A 10956. Left panel: Positive regions identified by thresholding, 
ranked in order of area. Center panel: Same for negative regions. Right panel: Corresponding 
STEREO /EUVI image. The smaller volume centered on the observed loop used for the 
energy calculations is illustrated by the box shown. 
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